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Abstract 

A convergence theorem for the continuous weak approximation of the solution of 
stochastic differential equations by general one step methods is proved, which is an 
extension of a theorem due to Milstein. As an application, uniform second order 
conditions for a class of continuous stochastic Runge-Kutta methods containing 
the continuous extension of the second order stochastic Runge-Kutta scheme due 
to Platen are derived. Further, some coefficients for optimal continuous schemes ap- 
plicable to ltd stochastic differential equations with respect to a multi-dimensional 
Wiener process are presented. 
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1 Introduction 



Since in recent years the application of stochastic differential equations (SDEs) 
has increased rapidly, there is now also an increasing demand for numerical 
methods. Many numerical schemes have been proposed in literature, see e.g., 
Kloeden and Platen jH] or Milstein and Tretyakov [11J and the references 
therein. The present paper deals with conditions on the local error of continu- 
ous one step methods for the approximation of the solution X(t) of stochastic 
differential equations such that global convergence in the weak sense is assured. 
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The continuous methods under consideration are also called dense output for- 
mulas [5]. Continuous methods are applied whenever the approximation Y(t) 
has to be determined at prescribed dense time points which would require the 
step size to be very small. Also if a graphical output of the approximate so- 
lution is needed, continuous methods may be applied as well. Further, for the 
weak approximation of the solution of stochastic differential delay equations 
with variable step size, a global approximation to the solution is needed. In the 
deterministic case it is known [5J that continuous methods may be superior 
over many other interpolation methods. Therefore, the application of continu- 
ous methods for stochastic differential delay equations may be very promising 
for future research (see also [T]). 

As an example, the continuous extension of a certain class of stochastic Runge- 
Kutta methods is given. Some stochastic Runge-Kutta (SRK) methods for 
strong approximation have been introduced by Burrage and Burrage j2~tf5] 
and Newton [12]. For the weak approximation, SRK methods have been de- 
veloped as well, see e.g., Komori and Mitsui [8J, Komori, Mitsui and Sug- 
iura [7], RoBler [T3fl4fl5|IT6] or Tocino and Ardanuy [T7] and Tocino and 
Vigo-Aguiar [T%] . 

The main advantage of continuous SRK methods is the cheap numerical ap- 
proximation of E(f(X(to + Oh))) for the whole integration interval < 6 < 1 
beside the approximation of E(f(X(to+h))). Here, cheap means without addi- 
tional evaluations of drift and diffusion and without the additional simulation 
of random variables. 

The paper is organized as follows: We first prove a convergence theorem for 
the uniform continuous weak approximation by one step methods in Section |2l 
which is a modified version of a theorem due to Milstein [TO]. Based on the 
convergence theorem, we extend a class of SRK methods of weak order one 
and two to continuous stochastic Runge-Kutta (CSRK) methods. Here, the 
considered class of SRK methods contains the well known weak order two SRK 
scheme due to Platen [6] . Further, we give order conditions for the coefficients 
of this class of CSRK methods in Section [3j Taking into account additional 
order conditions in order to minimize the truncation error of the approxima- 
tion yields optimal CSRK schemes. Finally, the performance of one optimal 
CSRK scheme is analyzed by a numerical example in Section HI 
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2 Local and Global Weak Convergence 



Let (X(t))tei denote the solution process of a d-dimensional Ito stochastic 
differential equation 

dX(t) = a(t,X(t))dt + b(t,X(t))dW(t), X(t )=X , (1) 

with a driving m- dimensional Wiener process (W(t)) t > and with I = [t ,T]. 

It is always assumed that the Borel-measurable coefficients 

and b : I x M d — y ]jj dxm satisfy a Lipschitz and a linear growth condition 

\\a(t,x) - a(t,y)\\ + || - b(t,y)\\ < C \\x - y\\, (2) 
||a(t,a;)|| 2 + ||6(t,x)|| 2 <C 2 (l + ||a;|| 2 ), (3) 

for every t G /, x, y G M. d and a constant C > such that the Existence and 
Uniqueness Theorem 4.5.3 [B] applies. Here, || • || denotes a vector norm, e.g., 
the Euclidian norm, or the corresponding matrix norm. 

Let a discretization Ih = {to, ti, • • • , ^tv} with to < t\ < . . . < t/v = T of the 
time interval J = [to, with step sizes h n = t n+ \ — t n for n = 0,l,...,iV — 1 
be given. Further, define h = max < n <jv h n as the maximum step size. Let 
X to ' x ° denote the solution of the stochastic differential equation (prj in or- 
der to emphasize the initial condition. If the coefficients of the Ito SDE (jTJ 
are globally Lipschitz continuous, then there always exists a version of X s,x (t) 
which is continuous in (s, x, t) such that X to > x °(t) = X s ' Xt °' x °^(t) P-a.s. holds 
for s,t G I with s < t (see [9]). For simplicity of notation, in this section it is 
supposed that Ih denotes an equidistant discretization, i.e. h = (% — t )/N. 
Next, we consider the one-step approximation 

Y t >*(t + h) = A(t,x,h;£), (4) 

where £ is a vector of random variables, with moments of sufficiently high 
order, and A is a vector function of dimension d. We write Y n+1 = Y tn,Yn (t n+ i) 
and we construct the sequence 

Yq = Xq, 

Y n+1 = A(t n , Y n , h;C n ), n = 0, 1, . . . , N - 1, U 

where £ is independent of Y , while £ n for n > 1 is independent of Y , . . . , Y n 
and £o,---,£n-i- Then Y = {Y s ' x (t) : x G R, s < t with s,t G 4} is a 
(non-homogeneous) discrete time Markov chain depending on h such that 
Y to ' Xo (t) = Y s > YkhX °( s \t) P-a.s. holds for s,t G 4 with s < t 

In the following, let Cj,(R d ,]R) denote the space of all g G C'(]R d ,IR) with 
polynomial growth, i.e. there exists a constant C > and r G N, such 
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that \d l x g(x)\ < C(l + ||x|| 2r ) holds for all x G R d and any partial deriva- 
tive of order i < I. Further, let g G C k P \l x R d , R) if g(-,x) G C k (I,R) and 
g(t, •) G Cp(R d , R) for all t G / and x G R d 0. 

Since we are interested in obtaining a continuous global approximation con- 
verging in the weak sense with some desired order p, we give an extension of the 
convergence theorem due to Milstein [TOffTT] which specifies the relationship 
between the local and the global approximation order. 

Theorem 2.1 Suppose the following conditions hold: 

(i) The coefficients a 1 and 6 M are continuous, satisfy a Lipschitz condition 
HP and belong to Cp +1,2(p+1) (/ x R d , R) with respect to x for i = 1, . . . ,d, 
j = l,...,m. 

(ii) For sufficiently large r (specified in the proof) the moments £'(||y n || 2r ) 
exist and are uniformly bounded with respect to N and n = 0, 1, . . . , N. 
(m) Assume that for all f G C 2 p {p+1) (R d , R) there exists a K G C P (R d ,R) 
such that the following local error estimations 

\E{f{X^ x {t + h)))- E{f{Y^ x (t + h)))\< K(x) h p+1 (6) 
\E(f{X^ x (t + 6h))) - E(f{Y^ x {t + eh)))\<K{x)h v (7) 

are valid for x G R d ; t,t + h G [t , T] and 6 G [0, 1] . 

Then for all t G [to, T] the following global error estimation 

\E(f(X to ' Xo (t))) - E(f(Y t0 ' x "(t)))\ <Ch p (8) 

holds for all f G Cp P+1) (R d , R), where C is a constant, i.e. the method 
has a uniform order of accuracy p in the sense of weak approximation. 

Remark 2.2 In contrast to the original theorem now the order of convergence 
specified in equation is not only valid in the discretization times t G Ih- 
Provided that the additional condition is fulfilled, the global order of con- 
vergence holds also uniformly for allt G [to,T]. 

Proof of Theorem 12.11 The proof extends the ideas of the proof of Theo- 
rem 9.1 in [10]. However, now all time points t G [to>^] have to be considered 
instead of only t G Ih like in [10] and an additional estimate is necessary in 
order to prove the uniform order of convergence on the whole time interval. 
Therefore, we consider the function 

u(s,x) = E(f(X(t k + 9h))\X(s) = x) 

for s,t k + 9h G J, x G R d and t k G Ih with s < t k and for 9 G [0, 1]. Due to 
condition @ u(s, •) has partial derivatives of order up to 2(p+ 1), inclusively, 
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and belongs to C P '(M. d , R) for each s G [io>£fc] [6lfT0] ■ Therefore, w satisfies 

X'-^t + /i))) - E(u(s, Y l - X {t + /i)))| < K u (x) h p+1 
uniformly w.r.t. s G [to,tfc] f° r some K u G Cp(M d ,M). Since Yq = Xq, we have 

E(f(x to > x °(t k + 0/1))) = emx*^'* ^ + eh))) - Eifix^ih + eh))) 

+ E(f(x t ^(t k + eh))). 

(9) 

Furthermore, since X ll > Yi {t k + Oh) = X t2 ' xtl ' Yl ^(t k + eh), we have 
E{f{X^{t k + Oh))) = E(f(X^ xH ' Yl ^\t k + 0/0)) - E(f(X^ Y %t k + Oh))) 

+ E(f(x t ^(t k + eh))). 

(10) 

By © and dTDJ we get 

E(f(X to ' Xo (t k + 0h))) 

=E(f(X^ xt0 ^\t k + Oh))) - E(f(X^(t k + Oh))) 
+ E(f(X^ xtl ^^\t k + Oh))) - E(f(X^(t k + Oh))) 

+ E(f(x t ^(t k + eh))). 

The procedure continues to obtain 

E(f(X to > Xo (t k + 0h))) 
fc— 1 

= x: (^(/(^• jr * i,n ^)(t fc + two)) - E(f(x^^(t k + 0/1)))) 

i=0 

+ E(/(x'- y H4 + ^))). 

Recalling that F'°' x °(t fc + Oh) = Y tk ' Yk (t k + Oh), this implies the identity 

E (f(X^ x °{t k + 0/0)) - E(/(F'«' x "(t fc + Oh))) 

fc— i 

= £ [E{E{f{X^^ Yi ^\t h + 0/0) | X'^+x))) 

i=0 

-S(S(/(X* i + 1 ' y * i ' ri ^+ 1 )(* Jfe + 0/0) I Y^fe+x)))) 
+ E(f{X*^ Y "(t k + 0/0)) - EifiY^fa + 0/0)). 

(11) 
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According to the definition of u(t,x), the Jensen inequality and (11 ip imply 

\E(f(X t0 > x °(t k + 9h))) - EUiY^itu + 0fc)))| 

fc— 1 

Y, (E(u(t w ,X t *> Y '(t i + fc))) - E{u{t l+l X^{U + fc))); 

i=0 

+ (EifiX^ih + 0fc))) - £7(/(Y'*> y *(* fc + 0fc)))) 
fe-i 

< £ E (\E(u{t i+1 , X^iU + fc)) - u{t i+1 , Y^(t t + fc)) | Y- 

i=0 ' 

+ e (^(/(x'^fo + ^)) - f(Y tk ' Yk (h + eh)) | n) |) . (12) 



Notice that the functions u(t, ■) and /, which belong to C f p P+1 ^(R d , R) and 
satisfy an inequality of the form and ([7P, also satisfy along with it a 
conditional version of such an inequality. Let K G Cp(IR d ,IR) such that for 
both u(t, ■) and / the inequalities (jSJ) and (J7j) hold. Thus there exist r G N 
and C > such that 

\K(x)\<C(l + \\x\\ 2r ) (13) 
holds for all x G R d . Then flTJ) together with (JT3]) imply for < 6>fc < fc the 
estimate 

|^(/(X*°^o( tfc + eh))) _ + ^)))| 

fc— l 

<^^(|^)^ +1 I) + ^(I^) Wl) 

i=0 
fc— 1 

< Y c(i + ^dl^H) h p+l + c(i + £(M 2r )) (#fc) p 

i=0 

<A;C(1+ max ^(ll^f 7 ')) h p+1 + C(l + £(||nf )) hP . (14) 

0<«<fc— 1 

Assuming that condition (Jul) holds precisely for this 2r and applying finally 
k = ^ in (ED, © is obtained. □ 

In the following, we assume that the coefficients a 1 and 5 M satisfy assump- 
tion ([!]) of Theorem 12.11 Further, assumption flu]) of Theorem 12.11 is always 
fulfilled for the class of stochastic Runge-Kutta methods considered in the 
present paper (see [T4"fTS] for details). 



3 Continuous Stochastic Runge— Kutta Methods 



As an example, we consider the continuous extension of the class of stochastic 
Runge-Kutta methods due to RoBler [T4"lll5|ll6] which contains the weak order 
two Runge-Kutta type scheme due to Platen [6]. The intention is to approx- 
imate the solution of the Ito SDE (TTJ on the whole interval / = [t ,T] in 
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the weak sense. Therefore, we define the d- dimensional approximation process 
Y by an explicit continuous stochastic Runge-Kutta method having s stages 
with initial value Y(t ) = X Q and 



Y(t n + 6h n )= Y(t n ) + a(t n + cPh n , HP) h T 

i=i 



s m 



i=l k=l 
s m 



+ E E ( h) + A (2 V) %0 b\t n + Ci (1) ^, # w 



(15) 



+ E E ( $ \o) h) + ft (4) (<0 '-$&) b k (t n + cfh n , m 
i=i fe,/=i 



for 6* G [0, 1] and n — 0, 1, . . . , N — 1 with stage values 



#f = K(i n ) + E ^ a(t n + cfh n , iff) h n 

3=1 
i—1 m 

3=1 r=l 

H\ k) = Y{t n ) + £ Ag } a(t n + cfh n , H$ 0) ) h n 

3=1 
3=1 

H\ k) = Y(t n ) + ^ 4 2) a(t n + cfh n , fff ) h n 

3=1 

+ E4 2) b k {t n + cfh n ,Hf ] )Jh n 

3=1 



for i = 1, . . . , s and k = 1, . . . ,m. Here, the weights a^, /?} G C([0, 1], R) are 
some continuous functions for 1 < z < s. We denote a(#) = (ai(9), . . . , a s (6>)) T , 
0(»O(0) = (/?j r) (fl), . . . , /?i r) (^)) T , G R s and A (<z) , G R sxs for < q < 2, 
1 < r < 4 and with A^f = B\f = for j > i, which are the vectors and ma- 
trices of coefficients of the SRK method. We choose = A^ q 'e for < q < 2 
with a vector e = (1, . . . , 1) T [H]. In the following, the product of vectors is de- 
fined component-wise. The coefficients of the CSRK method can be arranged 
in the following Butcher tableau: 
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c (0) 




5(0) 




c (l) 




5(1) 




c (2) 


^(2) 


£(2) 






T 

a 










/3( 3 ) r 


/3W T 



The random variables are defined by Im = AW 7 ^ and iWn = |(AW^ AW 7 ^ + 
V^''). Here, the All 7 ^ are independently three-point distributed random vari- 
ables with P(AW* = ±V3K) = | and P(AW£ = 0) = § for 1 < r < m. 
Further, the V£' are independently two-point distributed random variables 
with P(V k > 1 = ±h n ) = \ for jfe = 1, . . . , jfe - 1, V k ' k = -h n and V k > 1 = -V^ k 
for I = k + 1, . . . , m and k = 1, . . . , m [BJ. 

The algorithm works as follows: First, the random variables and 
have to be simulated for 1 < k, I < m w.r.t. the actual step size h n . Next, 
based on the approximation Y(t n ) and the random variables, the stage values 
H<-°\ and H {k) are calculated. Then we can determine the continuous 
approximation Y(t) for arbitrary t G [t n , t n+ i] by varying 9 from to 1 in for- 
mula (115p . Thus, only a very small additional computational effort is needed 
for the calculation of the values Y(t) with t G I\Ih- This is the main advantage 
in comparison to the application of an SRK method with very small step sizes. 



Using the multi-colored rooted tree analysis due to Rofiler [TITTS] . we de- 
rive conditions for the coefficients of the continuous SRK method assuring 
weak order one and two, respectively. As a result of this analysis, we obtain 
order conditions for the coefficients of the CSRK method (j!5)) which coincide 
for 6 = 1 with the conditions stated in [T5fT6] . The following theorem for 
continuous SRK methods is an extension of Theorem 5.1.1 in |16j . 



Theorem 3.1 Let a\b ij G C 3 /(I x 



for 1 < % < d, 1 < j < m. If the 



coefficients of the continuous stochastic Runge-Kutta method < T73)) fulfill the 
equations 



a(0f 







I. 

4. (/3«(0) T e) 2 = 
7. pV>{9) T BV>e 



2. f5 {4 \6fe 



5. ft 2) (6) 












3. ^\ef 

6. /?«(0) T 5«e = O 







for 9 G [0, 1] and the equations 



a{l) T A^\ 



10. (^\l) T e)(a T B^e) 



9. a(l) T (B^e) 2 = \ 
11. (/3 (1) (l) T e)(/3 (1) (l) T A«e) 
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12. 




(lfA^e = 


13. 


fr 2 \l) T B^e = l 


14. 




(lfB^e = 1 


15. 


(/3W(l) T e)(/3«(l) T (S«e) 2 ; 


16. 


{^\l) T e){^\l) T {B^e) 2 ) = ± 


17. 


/ g( 1 )(i) T (s( 1 )(s( 1 ) e )) = o 


18. 


/3< 3 > 


(l) T (B ( - 2 \B^e)) = 


19. 


£(3) (1) T (A (2) (jB (0) e))=0 


20. 




(l) T (A^(B^e)) = 


21. 


a(l) T (£(°)(£« e )) = 


22. 




(l) T A^e = 


23. 


/3( 4 )(l) T A( 2 )e = 


24. 


0(1) 


(l) T ((A«e)(B«e)) = 


25. 


0(3>(l) r ((A( 2 > e )(fl( 2 >e)) = O 


26. 


/3< 4 > 


(l) T (A^(B^e)) = 


27. 


^ 2 )(l) T (A«(5^e)) = 


28. 




(l) T (A^(B^e) 2 )=0 


29. 


^(l) T (A( 2 »(B (0) e) 2 ) = 


30. 




(l) T (B^(A^e)) = 


31. 


/3«(l) T (£«(A«e)) = 


32. 




(l) T (S«e) 2 = 


33. 


/^(l) T (£( 2 >e) 2 = 


34. 




(l) T (B ( - 2 \B^e)) = 


35. 


0( 2 >(l) r (BM(B« e )) = O 


36. 


0(1) 


(lf(B^e) 3 = 


37. 


/ g(3)(lf( S (2) e )3 = 


38. 


0(1) 


(l) T (5«(5«e) 2 ) = 


39. 


£(3) (1) T (jB (2) (jB (l) e)2)=0 



40. 


a(l) T ((B^e)(B^(B^e))) = 





41. 


^)(l) r ((iW(B(°) e ))(5We)) 




42. 


^(3) (1) T ((A (2 )(jB (o) e))(jB (2) e)) 




43. 


pm(l) T (AW(B(°\BWe))) = 





44. 


^(3) (1) T (A (2) (jB (0) (jB (l) e))) = 





45. 


/3«(1) T (5«(A«(5(%))) = 





46. 


0(3)(l) r (B( 2 )(A«(fl(°>e))) = 





47. 


/3( 1 )(l) r ((S( 1 )e)(B«(BWe))) 




48. 


/3 (3) (1) T ((jB (2) e)(jB (2) (jB (l) e))) 




49. 


/3«(l) T (5«(£M(5«e))) = 





50. 


^(3) (1) T (jB (2) (jB (l) (jB (l) e))) = 






are fulfilled then the continuous stochastic Runge-Kutta method < T73)) converges 
with order 2 in the weak sense. 



Proof. This follows directly from the order conditions for stochastic Runge- 
Kutta methods in Theorem 5.1.1 in (16] with the weights replaced by some 
continuous functions and taking into account 9 h instead of h for the expansion 
of the solution. Considering the order 1 conditions for the time discrete case, 
we obtain for example a T e h = h and ((3^ e) 2 h = h. Here, the left hand side 
results from the expansion of the approximation process while the right hand 
side comes from the expansion of the solution process at time t+h. Now, in the 
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continuous time case, in order to fulfill condition flTJ) of Theorem 12. II we replace 
the weights of the SRK method by some continuous functions depending on 
the parameter 9 and we consider the solution at time t + 9h. Thus, we obtain 
the order conditions a(9) T eh = 9h and {(3^ {9) T e) 2 h = 9 h. The remaining 
order 1 conditions can be obtained in the same manner. However, the order 
2 conditions need to be fulfilled only at time t + h, i.e., for 9 = 1, due to ([6]) 
of Theorem 12.11 Thus, we arrive directly at the conditions of Theorem 13.11 
Further, condition On} of Theorem 12. II is fulfilled for the approximations Y(t) 
at time t G Ih, see [15] . □ 

Remark 3.2 We have to solve 50 equations for m > 1 for schemes of order 
2. However in case of m = 1 the 50 conditions are reduced to 28 conditions 
(see, e.g. 115^16^ ). Further, explicit CSRK methods of order 2 need s > 3 
stages. This is due to the conditions 6. and 15., which can not be fulfilled for 
an explicit CSRK scheme with s < 2 stages. 

Remark 3.3 The conditions 1.-7. of Theorem VS. 1\ for 9 = 1 are exactly the 
order conditions for continuous SRK methods of weak order 1. 

As for time discrete SRK schemes, we distinguish between the stochastic and 
the deterministic order of convergence of CSRK schemes. Let Ps = P denote 
the order of convergence of the CSRK method if it is applied to an SDE and 
let pd with p D > p s denote the order of convergence of the CSRK method 
if it is applied to a deterministic ordinary differential equation (ODE), i.e., 
SDE ([TJ with 6 = 0. Then, we write (pd,Ps) for the orders of convergence in 
the following [IB] . 

Next, we want to calculate coefficients for the CSRK method (fT5|) which fulfill 
the order conditions of Theorem 13.11 Since the conditions of Theorem 13.11 
coincide for 9 = 1 with the order conditions of the underlying discrete time 
SRK method, we can extend any SRK scheme to a continuous SRK scheme. To 
this end, we need some weight functions depending on 9 G [0,1] which coincide 
for 9 = 1 with the weights of the discrete time SRK scheme as boundary 
conditions. Then, the conditions 1.-50. of Theorem 13.11 are fulfilled for 9 = 1. 
Further, the weight functions have to fulfill the conditions 1.-7. of Theorem 13. II 
also for all 9 G [0, 1]. For 9 = 0, the right hand side of ffT51) has to be equal 
to Y(t n ), i. e., the weight functions have to vanish for 9 = 0, which yields the 
second boundary conditions. As a result of this, if coefficients for a discrete 
time SRK scheme are already given, then the calculation of coefficients for 
a CSRK method reduces to the determination of suitable weight functions 
fulfilling some boundary conditions. 

As a first example, we consider the order (1, 1) SRK scheme with s = 1 
stage having the weights ct\ = = 1 and f3[ 2 ^ = (3[ 3 ^ = = 0. This is the 
well known Euler-Maruyama scheme. For the continuous extension, we have to 
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Table 1 

Coefficients of the optimal CSRK scheme CRDI1WM with p D = 2.0 and p s = 1.0. 

replace the weights a\ and fix by some weight functions aii, fil® G C([0, 1], R) 
for k = 1, . . . , 4. Due to Remark 13.31 we only have the boundary conditions 

ai (0) = P[ k) = 0, ai (l) = ^(1) = 1 and /3( 2) (1) = /^W = $>Q) = 
for k = 1, ... ,4, which have to be fulfilled for an order (1,1) CSRK scheme. 
Therefore, we can choose e*i(0) = ^{9) = 9 and p[ 2) (9) = (3[ 3) (9) = f3{ 4) {9) = 
for 9 G [0,1]. This results in the linearly interpolated Euler-Maruyama 
scheme, which is a CSRK scheme of uniform order (1, 1). 

Because there are still some degrees of freedom in choosing the continuous 
functions G C([0, 1],R) for k = 1,...,4, we can try to find some 

optimal functions in the sense that additional order conditions are satisfied. 
This results in CSRK schemes with usually smaller error constants in the 
truncation error. Therefore, these schemes are called optimal CSRK schemes 
in the following. For the already considered example of the order (1, 1) CSRK 
scheme with s = 1 stage, we can calculate an optimal scheme as follows: firstly, 
considering the continuous condition 1. of Theorem 13. 1[ we obtain for 9 G [0,1] 
the additional condition 

and from the continuous condition 4. that 

^\e) = Cl V9 (16) 

for some c\ G {—1, 1}. Further, the continuous conditions 2., 3. and 5. result 
in 

(3{ 2) (9) = 0, 0?\9) = O, /?i»=0, (17) 
for all 9 G [0, 1]. Now, all continuous weight functions are uniquely determined 
for the optimal order (1,1) CSRK scheme, which is still a continuous extension 
of the Euler-Maruyama scheme. 

Now, let us consider the order (2, 1) SRK scheme RDI1WM proposed in [4] 
with two stages for the deterministic part and one stage for the stochastic one. 
The coefficients of RDI1WM for A^°\ B®, and A^ 2 \ are given 

in Table Q] and the corresponding weights are 

fi = G,f], £ (1) = [i,o], 0® = F*> = pV> = [O,O]. 
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Table 2 

Coefficients of the optimal CSRK scheme CRDI2WM with p D = 2.0 and p s = 2.0. 

These coefficients are optimal in the sense that additionally some higher order 
conditions are fulfilled. For the continuous extension, we want to preserve the 
deterministic order 2 of the scheme. Therefore, beside the order one conditions 
1.-7. of Theorem 13.11 for 9 = 1 we have to take into account the continuous 
version of the deterministic order one condition which coincides with condition 
1. in Theorem 13. II Thus, we are looking for some weight functions aci, a 2 , /3[® G 
C([0, 1],R) fulfilling the boundary conditions c*i(0) = a 2 (0) = /3? } (0) = 0, 

<*i(l) = \, « 2 (1) = ^ /3! 1) (1) = 1, 

and p{ 2 \l) = /3f } (l) = /?J 4) (1) = for k = 1,...,4. Further, due to the 
continuous version of the deterministic order one condition, the equation 

ai(0)+a 2 (0) = 6 (18) 

has to be fulfilled for all 9 6 [0, 1]. In order to save computational effort we set 
Pi (9) = for k = 2, 3, 4 for all 9 G [0, 1]. Again, there are still some degrees of 
freedom in choosing the continuous functions ai,a 2 , (3{ 1] G C([0,1],R). Now, 
we can proceed exactly in the same way as for the order (1, 1) CSRK schemes. 
Considering additionally the continuous order conditions 2., 3., 4. and 5., we 
obtain again ( |T6l) and i fPTj) . Finally, we calculate from the continuous version 
of condition 8. that 

a 2 (9) = ^ 2 

and thus c*i(0) = 9 -\9 2 due to (JTHD for all 9 G [0, 1]. The coefficients of the 
CSRK scheme CRDI1WM are given in Table EJ 

We consider now the continuous extension of the order (2, 2) SRK scheme 
RDI2WM with s = 3 stages considered in [4]. The coefficients of the SRK 
scheme RDI2WM for A®\ B®, an d , B^ can be found in 
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Table |2] and the weights are given as 



^ [2' 2'^]' [4' 8' 8-1 ' ^ Pi 4 i 4]' 

/3 (3) = [-4-,U]> ^ = [0,%-%. 

We proceed in the same way as for the scheme CRDI1WM by taking into 
account some additional order conditions. As the stage number of the deter- 
ministic part of the scheme is only two, we set 0.3(6) = 0. From the order 
1 conditions of Theorem O follows for a h /?f } G C([0,1],R), 1 < i < 3, 
1 < k < 4, that 

a x (0) = 6 - a 2 {6), 

$\o) = Cl Ve-2$\e), &\e) = $\o), 

(6) = -& 2) (9)- fli 2) (6), (3{ i} [6) = (6)- /3| 4) (6). 

(k) 

Further, the functions at and p\ have to fulfill for 6 = and 6 = 1 the 
boundary conditions 

*(i) = 5, A W (D = |, /f>(D = ^, /f>(D = -^, 
= i ^(i) = j£ /#>(!) = 

and «i(0) = (0) = for 1 < i < 3 and 1 < k < 4. Here, the continuous 
versions of the order conditions 8.-50. of Theorem 13.11 i.e., conditions 8.-50. 
with a(6) and /3 (fe) (6>) instead of a(l) and f3 (k \l) for k = 1, . . . , 4, and with 
appropriate right hand sides are automatically fulfilled except of conditions 
8., 9., 10., 11., 13., 14., 15., 16., 22., 32. and 33. The continuous versions of 
condition 8. and 9. yield 

a 2 (6) = \0\ (19) 
However, condition 10. alternatively yields 

a 2 (6) = \oK 

so one has to decide whether 8. and 9. or 10. should be fulfilled. The conditions 
22. and 32. coincide and provide 

$\6) = -&\d). (20) 

Further, condition 33. results in 

4 4) (0) = -M>). (21) 
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Table 3 

Coefficients of the optimal CSRK scheme CRDI3WM with p£> = 3 and ps = 2. 
Taking into account condition 13., we calculate that 



?( 2 )/ 



(22) 



and due to condition 14. that 



1 



(23) 



Now, if we combine conditions 32. and 33. with 13. and 14. then we can 
determine and uniquely as 



1 



0. 



2\/2 ' 2v/2 
Considering condition 11., which coincides with 15., one obtains 



(24) 



(25) 



Finally, condition 16. yields that 



(26) 



Thus, one can choose from these additional conditions in order to minimize 
the truncation constant. Especially, if we consider equations ( IT9|) - (l26|) then 
we obtain the scheme CRDI2WM presented in Table [2J 

If we allow now three stages for the deterministic part, i.e., 0:3 7^ 0, then we 
can construct explicit CSRK schemes of order p D = 3 and ps = 2. Therefore, 
we extend the order (3, 2) explicit SRK scheme RDI3WM introduced in [4] 
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with coefficients A®, B^°\ A®, and A^ , 5( 2 ) from Table [3] and with 
the weights 

« = [2 I 41 5(1) _ rl 3 31 5(2) = r Q V6 _V61 

Lg> 3' gJ' ^ U' 8' 8J ' " \- 1 4 ' 4 J ' 

/3 (3) = [-iUL /3 (4) = [0,^,-^]. 
Here, the well known order three conditions 

a T (A^e) 2 = \, a T (A^(A^e)) = -. 

3 O 

for deterministic ODEs are fulfilled [5]. Obviously, the scheme RDI3WM differs 
from RDI2WM only in A^°\ and a. Therefore, we choose the functions 
f3\ k) G C([0,1],R) for 1 < i < 3 and 1 < k < 4 in the same way as for 
the CSRK scheme CRDI2WM. Thus, we only have to specify the functions 
aij G C([0, for 1 < i < 3. For order = 3, the continuous version of 
condition 8. in Theorem 13.11 

a(9) T A^e = ~0 2 
has to be fulfilled as well. As a result of this, we yield the conditions 
a x {6) = ~a 3 (0) + 6-6 2 , a 2 (9) = 9 2 - ~ot z {0), 



with boundary conditions 03(0) = and 0:3(1] 



9 • 



In order to specify a 3 G C([0, 1],R), we consider the continuous version of 
condition 9. which has the unique solution 

a 3 (0) = a 3 (l)6 2 . (27) 

Alternatively, we can consider the continuous version of condition 10. which 
yields 

2 2 (2\/l5 - 9) + 70§ 



o 3 (0) 



9 v/15 -1 



However, instead of considering condition 9. or 10., we can also take into 
account the continuous version of the deterministic order 3 condition 

a{6) T {A^\A^e)) = - 

6 

which yields 

a 3 (6) = a 3 (l)9 3 . 

On the other hand, if we consider the continuous version of the deterministic 
order 3 condition 

a(6) T (A^e) 2 = \e 3 
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Table 4 

Coefficients of the optimal CSRK scheme CRDI4WM with po = 3 and p$ = 2. 
then we get for the deterministic part of the continuous SRK scheme that 

a*(0) = fp-±P. (28) 

So, one can choose from these additional conditions in order to minimize the 
truncation error of the considered scheme. For the continuous extension of the 
SRK scheme RDI3WM, we choose condition (12 7p which yields the coefficients 
of the CSRK scheme CRDI3WM presented in Table EJ 

Analogously to the procedure for the CSRK scheme CRDI3WM, we extend 
the order (3,2) explicit SRK scheme RDI4WM @] with coefficients A (0) , B {0 \ 
v4 (1) , and A (2 \ from Table Hand with the weights 

a = a 2 II m = [1 1 31 5(2) = r Q VE _V61 

1-6' 3' 6-1' r 1 L4' 8' 8-1 ' " L w ' 4 ' 4 J ' 

/3 (3) = hiil], /3 (4) = [0,^,-^]. 
Then, we obtain the conditions 

a x {9) = a 3 (9) + 9-9 2 , a 2 (6) = 9 2 - 2a 3 (9), 

with boundary conditions 0:3(0) = and 03(1) = |. Choosing again condition 
(1271) yields the coefficients of the CSRK scheme CRDI4WM given in Table H 

However, if we consider condition (12 8p instead of (|27|) . then we obtain co- 
efficients for the CSRK scheme CRDI5WM. The coefficients of the scheme 
CRDI5WM coincide with the coefficients of the scheme CRDI4WM, with the 
exception of the weights oti{9). From condition (|28|) . we calculate the weights 
ai (9) = p 3 - \9 2 + 9, a 2 {9) = 29 2 - ffl 3 and a 3 (9) = §# 3 - \9 2 for the scheme 
CRDI5WM. 
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4 Numerical example 



Now, the proposed CSRK scheme CRDI3WM is applied in order to analyze 
the empirical order of convergence. The first considered test equation is a 
linear SDE (see, e.g., [6]) 

dX(t) = aX(t)dt + bX(t)dW(t), X(0) = x , (29) 

with a,b,Xo £ R. In the following, we choose f(x) = x and we consider the 
interval I = [0, 2]. Then the expectation of the solution at time t £ / is given 
by E(X(t)) = x ■ exp(ai). We consider the case a = 1.5 with b = 0.1 and 
x = 0.1. 

As a second example, a multi-dimensional SDE with a 2-dimensional driving 
Wiener process is considered: 

+ \^X\t) ±X\t) + ±X\t)) {w\t)) 

with initial value X(0) = (1, 1) T and I = [0,4]. This SDE system is of spe- 
cial interest due to the fact that it has non-commutative noise. Here, we are 
interested in the second moments which depend on both, the drift and the 
diffusion function (see [6] for details). Therefore, we choose f(x) = (x 1 ) 2 and 
obtain 

E(f(X(t))) = exp(-t) . 

For t 6 /, we approximate the functional u = E(f(Y(t))) by Monte Carlo 
simulation using the sample average um,h = j|Em=i/( f '' ra '(i)) of indepen- 
dent simulated realizations Y^ m \ m = 1,...,M, of the considered approxi- 
mation Y and we choose M = 10 9 . Then, the mean error is given as fi = 
UM,h — E(f(X(t))) and the estimation for the variance of the mean error is 
denoted by a 2 . Further, a confidence interval [fi — A/t, fi + A/t] to the level of 
confidence 90% for the mean error fi is calculated [B]. 

First, the solution E(f(X(t))) is considered as a mapping from / to R with 
t i — y E(f(X(t))). Here, the whole trajectory of the expectation even between 
the discretization points has to be determined. Therefore, we apply the CSRK 
scheme CRDI3WM with step size h = 0.25 and determine E(f(Y(t n ))) for the 
discretization times t n = n ■ h for n = 0, 1, . . . , iV and the approximation Y(t) 
is exploited between each pair of discretization points t n and t n+ i by choosing 
9 £ ]0, 1 [ . This has been done for 9 = 0.1, . . . , 0.9. The results are plotted in 
the left hand side of Figure [T] and Figure |2J Further, the errors are plotted 
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t Id h 



Fig. 1. Scheme CRDI3WM with SDE ([29}. 




t Id h 



Fig. 2. Scheme CRDI3WM with SDE (|3TJjh 

along the whole time interval / in the figure below. Here, it turns out that the 
continuous extension of the SRK scheme works very well. 

Next, SDE ( 1291) and SDE ( |30l) are applied for the investigation of the or- 
der of convergence. Therefore, the trajectories are simulated with step sizes 
2-\ . . . , 2" 4 for SDE flU and with step sizes 2 1 , . . . , 2" 2 for SDE (j50|>. As an 
example, we consider the error fi at time t = 1.7 for SDE (1291) and at t = 3.8 
for SDE ( J30l . which are not discretization points. The results are plotted on 
the right hand side of Figure Q] and Figure [2] with double logarithmic scale 
w.r.t. base two. On the axis of abscissae, the step sizes are plotted against the 
errors on the axis of ordinates. Consequently one obtains the empirical order 
of convergence as the slope of the printed lines. In the case of SDE (T29]) we 
get the order p w 2.67 and in the case of SDE (1301) we get the order p w 2.27. 
Table and Table contain the corresponding values of the errors, the vari- 
ances and the confidence intervals to the level 90%. The same has been done 
for all time points in order to determine the empirical order of convergence 
on the whole time interval /. These results are listed in Table [7] for both con- 
sidered examples. The very good empirical orders of convergence confirm our 
theoretical results for the CSRK scheme CRDI3WM. 
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Table 5 

Scheme CRDI3WM with SDE 



23) at i = 1.7. 



Table 6 
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A 




jl — Ajl 


jl + Ajl 




-2.188E-02 


1.146E-09 


-2.189E-02 


-2.187E-02 


2- 2 


-3.965E-03 


1.705E-09 


-3.975E-03 


-3.956E-03 


2- 3 


-5.662E-04 


1.715E-09 


-5.760E-04 


-5.564E-04 


2 _4 


-8.682E-05 


1.746E-09 


-9.672E-05 


-7.691E-05 


DI3WM with SDE ( 


3D) at t = 3.8. 


h 


A 




jl — Aft 


ji + Ajl 


2 1 


-1.031E-02 


9.365E-12 


-1.031E-02 


-1.031E-02 


2° 


-2.161E-03 


2.294E-11 


-2.162E-03 


-2.160E-03 


2- 1 


-4.258E-04 


4.453E-11 


-4.274E-04 


-4.242E-04 


2- 2 


-9.392E-05 


3.375E-11 


-9.530E-05 


-9.255E-05 



Table 7 

Orders of convergence for CRDI3WM. 
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SDE (|30D 
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order 
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order 
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order 
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order 


0.1 


2.35 


1.1 


2.94 




0.2 


1.82 


2.2 


2.30 


0.2 


3.39 


1.2 


2.60 




0.4 


1.84 


2.4 


2.27 


0.3 


3.54 


1.3 


2.69 




0.6 


1.95 


2.6 


2.24 


0.4 


3.16 


1.4 


2.67 




0.8 


2.07 


2.8 


2.24 


0.5 


2.87 


1.5 


2.80 




1.0 


2.13 


3.0 


2.22 


0.6 


3.50 


1.6 


2.89 




1.2 


2.26 


3.2 


2.25 


0.7 


2.32 


1.7 


2.67 




1.4 


2.32 


3.4 


2.26 


0.8 


2.44 


1.8 


2.72 




1.6 


2.36 


3.6 


2.27 


0.9 


2.53 


1.9 


2.70 




1.8 


2.37 


3.8 


2.27 


1.0 


2.79 


2.0 


2.80 




2.0 


2.31 


4.0 


2.24 
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